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Abstract 

Sizing methods can be used to get a first estimate of the required Hall thruster dimen- 
sions and operating conditions for a given input power and a corresponding thrust level. 
After a review of the existing methods, a new approach, which considers the three char- 
acteristic thruster dimensions, i.e. the channel length, the channel width and the channel 
mean diameter as well as the magnetic field, is introduced. This approach is based on 
analytical laws deduced from the physical principles that govern the properties of a Hall 
effect thruster, relying on a list of simplifying assumptions. In addition, constraints on the 
channel wall temperature as well as on the propellant atom density inside the channel are 
taken into account. The validity of the scaling laws is discussed in light of a vast database 
that comprises 23 single-stage Hall effect thrusters covering a power range from 10 W to 
50 kW. Finally, the sizing method is employed to obtain a preliminary geometry and the 
magnetic field strength for a 20 kW and a 25 kW Hall effect thruster able to deliver a 
thrust of IN, respectively 1.5 N. 

To be published in Journal of Technical Physics 49, vol. 3-4 (2008). 
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Nomenclature 



A channel cross section; area 

B magnetic field strength 

C/^p analytic proportionality coefficient for I^p 

Cj^^ empiric proportionality coefficient for Igp 

Cti 2 analytic proportionality coefficient for T 

C^^ 2 empiric proportionality coefficient for T 

Cp analytic proportionality coefficient for P 

Cp empiric proportionality coefficient for P 

d, dext, dint mean, external and internal channel diameter 

Gij grey body configuration factor 

h channel width 

Id, li discharge current, ion current 

lyji current corresponding to the propellant mass flow rate 

Igp specific impulse 

kB Boltzmann constant 

L channel length 

mn,mi,me propellant atom mass, ion mass, electron mass 

m, rhi propellant mass flow rate through the anode, ion mass flow rate 

nn,ne,ni atom, electron and ion number density 

P input power 

Pioss power losses due to plasma-wall interactions 

Qp heat flux deposited by the plasma onto channel walls 

r, R inner, respectively outer, channel radius 

fLe^fLi electron and ion Larmor radius 

T thrust 

Tmax, Text, Tint maximum, external and internal wall temperature 

Tf, electron temperature 

Ud discharge voltage 

Vn,Ve,Vi thermal velocity of atoms, electrons and ions 

Vi average ion exhaust velocity 

a propellant conversion efficiency 

P ratio between ionization mean free path and channel length 

A voltage losses 

e emissivity 

r] thrust efficiency 

Aj , Aee ionization and electron-electron collision mean free path 

Vce electron gyrofrequency 

electron-atom collision frequency 

^ scaling index variable 

(Ji,aee cross-section for ionization and electron-electron impact 

a^n electron-atom momentum exchange cross-section 

Tc gyroperiod 

Ten electron-atom coUisional time 
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1 Introduction 



Electric propulsion is nowadays a well-established technology for space applications pP . Among 
all proposed electric propulsive devices such as arcjet, magnetoplasmadynamic thruster, grid- 
ded ion engine and Hall Effect Thruster (HET), the latter is currently recognized as an 
attractive propulsion means for long duration missions and for maneuvers that require a large 
velocity increment. Hall effect thrusters, also called Stationary Plasma Thrusters or closed 
electron drift thrusters, are advanced propulsion devices that use an electric discharge with 
magnetized electrons to ionize and accelerate a propellant gas [21 [3] . Due to interesting fea- 
tures in terms of propellant ejection speed, efficiency, flexibility and lifetime, HETs are now 
employed for missions like geo-stationary satellite orbit correction and station keeping. More- 
over, HETs appear as good candidates to be used as the primary propulsion engine for space 
probes during interplanetary journeys, as demonstrated by the successful SMART-1 Moon 
flyby solar-powered mission of the European Space Agency [1]. 

A schematic of a Hall effect thruster is depicted in Fig. [TJ The basic physics of a HET 
implies a magnetic barrier and a low pressure DC discharge generated between an external 
hollow cathode and an anode in such a way that a crossed electric and magnetic fields discharge 
is created [UEIE]. The anode, which also serves as gas injector, is located at the upstream 
end of a coaxial annular dielectric channel that confines the discharge. Xenon is generally 
used as a propellant gas for its specific properties in terms of high atomic mass and low 
ionization energy. A set of solenoids provides a radially directed magnetic field B of which 
the strength is maximum in the vicinity of the channel exhaust. The magnetic field is chosen 
strong enough to make the electron Larmor radius much smaller than the discharge chamber 
length, but weak enough not to affect ion trajectories. The electric potential drop is mostly 
concentrated in the final section of the channel owing to the low electron axial mobility in this 
restricted area. The corresponding induced local axial electric field E has two main effects. 
First, it drives a high electron azimuthal drift - the Hall current - that is responsible for the 
efficient ionization of the supplied gas. Second, it accelerates ions out of the channel, which 
generates thrust. The ion beam is neutralized by a fraction of electrons emitted from the 
hollow cathode. When operating near 1.5 kW, a HET ejects ions at 20kms~^ and generates 
100 mN of thrust with an overall efficiency of about 50 %. 

New fields of application are nowadays envisaged for electric propulsion systems that 
require low and high power devices. Low power Hall thrusters (~ 100 W) are well suited 
for drag compensation of observation satellites that operate on a low-altitude orbit in the 
Earth atmosphere as well as for orbit correction and orbit transfer of small platforms. The 
use of high power (~ 5 kW) Hall thrusters for orbit raising and orbit topping maneuvers of 
communication satellites would offer significant benefits in terms of launch mass, payload mass 
and operational life. In addition, journeys towards far-off planets and asteroids with large 
and heavy robotic probes necessitate to build thrusters with an input power in the range 
10-100 kW. In view of the projects demand, it appears necessary to expand the operating 
envelope of existing Hall effect thruster technology to achieve the required performance level. 
A non-trivial question then arises: How to extrapolate the design and architecture of currently 
existing Hall thrusters towards different scales and different operating conditions? In other 
words, what are the scaling laws that connect Hall effect thruster characteristic dimensions 
with operating parameters like discharge voltage, propellant mass flow rate and magnetic field 
strength and performances in terms of thrust, specific impulse and overall efficiency? 

Scaling laws that govern the physical properties, the accelerating ability as well as propel- 
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lant and energy consumption of Hall thrusters have been extensively investigated by numerous 
authors since the period of development of Hall thrusters in the 70's. In spite of decades of 
research on this subject, the assessment of scaling laws is still a topic of great interest with 
debates and controversies as various methodologies and results exist. Therefore, before de- 
scribing the approach associated with this study as well as its outcomes, it is worth briefly 
reviewing works carried out on this topic during the past few years and available in the lit- 
erature. According to the pioneer works of the Russian physicist A. Morozov [5], in order 
to derive scaling laws, it is necessary to find a similarity criterion, or a set of criteria, that 
warrant the underlying physical processes stay unchanged whatever the truster. This prin- 
ciple states that the properties of thrusters with a different geometry are linked by way of 
scaling laws only if a sufficiently large number of dimensionless similarity criteria coincide. 
The complete set of similarity criteria has not yet been found, however, all works show that 
the Melikov-Morozov criterion has a strong impact on HET behavior and characteristics and 
it must always be taken into account. This criterion indicates that A,/ <^ L. In addition to 
similarity criteria, the investigation of scaling laws for Hall thrusters necessitates to account 
for simplifying assumptions, e.g. a frozen electron temperature, and constraints like a high 
efficiency and a reasonable channel wall temperature. 

V. Zhurin and co-workers propose a size scaling method which is limited to the effects of 
changing either the channel width h or the channel mean diameter d [2]. One assumption 
is that the channel mean diameter is much larger than the channel width, so that variation 
of parameters in the radial direction are negligible. Furthermore, they considered a fixed 
discharge voltage. In order to obtain similar performances for two different thruster con- 
figurations, authors shown that the ratio of r^, A and L to h must stay the same for the 
two configurations. Using these criteria and a set of assumptions they demonstrate that the 
magnetic field strength is inversely proportional to the channel width, B 1 /h, whereas the 
discharge current and the propellant mass fiow rate are inversely proportional to the channel 
mean diameter, oc 1/d and m cx 1/d. 

The method presented by J. Ashkenazy et al. concerns low power scaling and it is based 
on the idea of a sufficient propellant utilization [7]. They show by means of a simplified 
analysis that a straightforward approach for scaling down the channel size results in a rise of 
power losses and a reduced overall thruster efficiency. To avoid these effects, the acceleration 
region has to be scaled down along with the channel width and the magnetic field strength 
must be increased in proportion to 1/h. They also describe an alternative view that consists 
in extending the channel length. 

M. Martinez-Sanchez et al. propose an approach for low-power Hall effect thrusters that 
includes the use of reference points j^. The goal is to achieve a reduction in the thrusters 
length scale while preserving both the thrust efficiency and the specific impulse. The main 
result of their study is that the propellant mass flow rate and the applied power scale as 
the channel length, m oc L and P oc L, whereas the magnetic field strength is in inverse 
proportion to the latter, B oc 1/L. This scaling method allows to calculate the channel size 
and the performances of a small thruster with a given input power provided that a well-known 
thruster is utilized as a reference. 

The team led by M. Andrenucci suggests to employ an analytical model coupled to an 
existing Hall thruster database [9]. This scaling methodology aims to provide design options 
for high power thrusters. The idea is to create a vector of fundamental parameters defining 
the thruster geometry and its performances. This vector is composed of three geometric 
parameters (L, h and d), the gas number density in the injection plane and the applied 
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discharge voltage Ud- A scaling matrix derived from the Hall thruster physical principles 
allows to obtain a new thruster characteristics on the basis of a reference thruster. 

In a recent article [10], Daren Yu et al. present an improvement of the existing scaling 
theory by introducing a scaling index variable. They assume that the electron temperature 
and the discharge voltage are constant, that the ratio Xi/L is constant and that the geometric 
similarity is given by r/i?=constant and i?^^^/L=constant, where R and r are the channel 
outer respectively inner radius and ^ represents a scaling index variable. They show by way 
of a comparison of experimental data with numerical outcomes for different values of ^, that 
results obtained from their scaling theory agree well with the experimental data for .^=2. 
Therefore they deduce that the number density n is constant, whereas the mass flow rate rh, 
the input power P, the thrust T and the discharge current Id are proportional to the square 
of the outer channel radius R. 

Finally, it was demonstrated by V. Kim [B] that to reach a high level of efficiency it is 
not only necessary to ionize and accelerate ions but to accelerate them into the proper direc- 
tion, hence the need for an optimized magnetic field topology. In short, for a HET with an 
optimized magnetic field map, there is a relationship between the acceleration layer length 
La and the magnetic layer characteristic size Lb and between La and h. The use of these 
similarity criteria, which include the magnetic field structure, permits to simplify the devel- 
opment of new Hall thrusters models. V. Kim et al also emphasize the fact that in the case 
of a traditional HET design, the ratio of the ionization mean free path to the channel length 
must stay constant [Melikov-Morozov criterion) [11]. Moreover they give additional criteria 
about the geometry of thruster elements, i.e. L and d are both proportional to /i, that are 
equivalent to the ones given by V. Zhurin. There are two immediate consequences. First, 
the propellant mass flow density and the power density rise when the thruster size decreases, 
assuming a constant discharge voltage. Second, as the magnetic field strength is connected 
with the characteristic dimensions, notably the channel width, and with the operation mode, 
the magnetic field magnitude must rise when the size reduces in order to keep comparable 
conditions. 

In this contribution, we propose an original way to extrapolate Hall thruster geometry 
towards both the low and high power ranges. The approach is based on the combination of 
a set of scaling laws with a vast body of existing data. Besides, realistic constraints on the 
performance level and on the thermal load allow to refine the set of possible L, h and d char- 
acteristic scale lengths found for a given thruster input power. In section 2, similarity criteria 
and scaling laws are derived from the physics at work in a HET and, consequently, rela- 
tionships between performances, operating parameters and the characteristic dimensions are 
established. The thermal constraint is discussed in section 3 along with modeling of thruster 
thermal behavior. A detailed description of the database on Hall thruster performances and 
properties is given in section 4. Besides, recorded data are compared with calculation out- 
comes. In section 5, our approach is finally employed to design a 20 kW as well as a 25 kW 
thruster. A summary and general conclusions follow in section 6. 



5 



2 Set of governing relations 



2.1 Similarity criteria and scaling laws 

A necessary first step in order to determine scaling laws for Hall effect thruster does consist in 
finding some critical parameters as well as in defining similarity criteria based on the current 
knowledge and understanding of the physics of Hall thrusters. The geometry of a Hall effect 
thruster is defined by three characteristic dimensions, see Fig. [H the discharge channel length 
L, the mean channel diameter d = \{dext + c?mt) and the channel width h, as well as by a set 
of operating parameters such as the magnetic field strength B, the discharge voltage and 
the propellant mass flow rate rn. 

To simplify the assessment of scaling laws, the following assumptions have been made through- 
out the entire paper: 

- the electron temperature is constant and homogeneous whatever the operating condi- 
tions, 

- the propellant gas has a uniform temperature all over the channel, hence a constant 
propellant velocity, 

- the potential energy is fully converted into kinetic energy and all ions experience the 
whole potential drop, of which the magnitude is C/^, 

- plasma-wall interactions are taken into account through heat load to the channel walls, 

- the magnetic field is uniform; its value at the exit plane is solely considered, 

- electron transport across the magnetic barrier is considered as classical: no anomalous 
transport is accounted for within the region of strong magnetic field [12j . 

- there are no multiply-charged ions in the plasma, 

- a parallel monokinetic ion beam is produced, i.e. the plasma jet divergence is null. 

The first relationship reflecting the impact of the thruster scale on its performance is the 
relationship between the discharge channel length L and the ionization mean free path Aj. 
To ensure a sufficient ionization of the gas, it is necessary to satisfy the Melikov-Morozov 
criterion: 



(2.1) 



Xi <. L. 



The ionization frequency that originates from electron-atom impacts reads: 



(2.2) 




The ionization length, which corresponds to the mean distance an atom can travel before being 
ionized, can be formulated as the product Ui u„ in a first order approximation. Therefore the 
Melikov-Morozov criterion can be expressed as: 



Assuming that both the electron temperature Tg and the propellant speed Vn are constant 
and independent of the operation conditions, imphes that ai and Ve, being only a function of 
Te, remain constant. To keep the ratio Xi/L constant, equation 12.31 leads to: 

(2.4) rin^j. 

In the same way, the ratio Xee/L must stay constant, where Age is electron-electron impact 
mean free path: 

(2.5) Aee = , 

where the cross-section aee is a function of the Coulomb logarithm [H] . The equation 12.51 
implies that: 

(2.6) Ue oc — . 

The relations 12.41 and 12.61 for the atom, respectively the electron, number density are the 
same than those developed before by V. Zhurin and M. Martinez- Sanchez [21 [8]. Another 
relation between n„ and Hall thruster dimensions can be established when considering the 
propellant mass flow rate passing through the anode. The quantity m can be decomposed 
into the product of several terms: 

(2.7) rh = rin ■ rrin ■ Vn ■ A. 
The annular channel cross section A is given by: 

(2.8) A = y(c?ext - dint) = T (dext + dint) (dext " dint) = T^dh, 

4 4 V ^ 'V ^ ' 

2d 2h 

hence A oc hd. We can therefore consider that for a constant gas mass flow rate: 

(2.9) un oc ^. 

This second relationship between the atom number density and the thruster dimensions h and 
d has never been mentioned previously, as authors never considered two sizes at the same time. 
Note that, in order to keep the physical processes at work in a Hall effect thruster unchanged, 
to warrant a high efficiency and to limit the thermal load as well as the channel wall wear, the 
number densities of electron and atoms must stay roughly fixed inside the thruster channel 
whatever the input power. The average values commonly found in literature, which turn out 
to guarantee a satisfying operation, are: ?i„ = 10^^ m~^ and Ue = 10^^ m~^. 
Another interesting fact is that the variables L, h an d are linked to one another. Indeed, by 
using the Melikov-Morozov criterion and the mass flow rate definition, it comes: 

p =^ = pL =^ n. 



L nnCTiVe L aiVeP 

Lj ml 
m = n„ ■ m„ ■ v„ ■ nnd iir, 



m„f„7r hd 
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Combining the two previous relations, one finds: 



2.10 — = 2~ = /(/^' 



Therefore, for a fixed value of the [3 parameter, and under our list of assumptions, the thruster 
characteristic dimensions are coupled through the gas mass flow rate. 

The magnetic field strength in a Hall effect thruster is such that electrons are magnetized and 
ions are not magnetized. The following criterion must then be fulfilled [6]: 



(2.11) rie^L^ rii 
The definition of the electron Larmor radius is: 

(2.12) TLe 



eB 

Using the fact that the ratio of rie to L must remain constant, the following relationship can 
be established: 

(2.13) Socy. 

The relation 12.131 between B and L has already been mentioned by M. Martinez-Sanchez [8]. 
A second constraint for the magnetic field strength can be established due to the fact that the 
electron gyroperiod Tee in the magnetic barrier must be much shorter than the time between 
two consecutive electron-atom collisions: 

(2.14) !£!! = » 1, 

Tee ^en 

where fee = eB/nie and u^n = nnC^enVe [IS]- This strong point indicates that electrons must 
be efficiently trapped inside the magnetic field of a Hall thruster in order to produce a high 
electric field and to favor ionization of the seeded gas.. In fact, Ten is so long in a HET that 
anomalous electron transport perpendicular to the magnetic field lines must be put forwards 
to correctly explain the observed properties and the magnitude of measured quantities [21 [5] . 
Assuming that the ratio Vcl^en must remain constant. Eg. 12.14] implies: 

(2.15) ^'^Yd 

To the best of our knowledge, the correlation between B and a product of two dimensions 
has never been mentioned before as solely one characteristic thruster dimension is usually 
considered by authors. Combining equations 2.9 and 2.15, one obtains i? tx n„ in compliance 
with the fact that plasma containment depends on collision events with neutrals. 



2.2 Relationship between performances and dimensions 

The previously derived scaling laws can now be utilized to analyze the effect of dimensions 
upon the performances of a Hall effect thruster. 
The specific impulse is defined as follows: 

(2.16) /sp = -, 

50 
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where go is the standard gravity at Earth's surface. The average ion exhaust velocity can be 
defined as [3]: 



(2.17) Vi = J—{Ud-A). 

V m 

Voltage losses are neglected in this work. Therefore the ion exhaust velocity is proportional 
to the square root of the discharge voltage {vi tx and for this reason: 



(2.18) Isp a VUd- 
The thrust of a Hall effect thruster reads: 

(2.19) T = rriiVi = arhvi, 

where the coefficient a stands for the fraction of propellant atoms that are converted into ions. 
Considering in a first approach that a is constant - the typical value given in the literature is 
a ~ 0.9 - and using equation (j2.7p . one can establish the following relationship for the thrust: 

(2.20) T (X y y^hd. 

Assuming that the discharge current has no electronic component, i.e. h, one can write: 

(2.21) Ifi a ■ m cm a ■ —hd. 

ij 

These relations between the performances and the dimensions are in agreement with those 
described by Daren et al. [10], if one considers, like they do, a constant number density n„, 
a constant discharge voltage Ud and a geometric similarity R (xr. Using the relation 12.211 for 
the discharge current, one finds out that the applied power P = Udid depends on the thruster 
characteristic dimensions. In like manner, the thrust efficiency is expected to depend on the 
thruster geometry. The efficiency is defined as: 

(2.22) r] 



2mUdId 



The efficiency must account for the gas flow rate injected through the anode and supplied to 
the cathode. Most of the time only the gas flow injected into the channel is considered; one 
then talks about the anode efficiency. Surprisingly, due to the chosen assumptions, the anode 
thrust efficiency r/ is solely a function of the propellant conversion coefficient a. 



3 Thermal constraint 

During thruster operation, a certain percentage of the input power P is lost due to plasma- 
wall interactions. Indeed, as shown in [16], a relatively large energy flux qp is deposited onto 
the discharge channel walls, mostly due to ion and electron bombardment, which results in 
a temperature increase of all thrusters components. Naturally, there is maximum amount of 
power that can be passed to the walls in order to limit the thermal load and to minimize the 
sputtering yield of the wall material, usually BN-Si02. On can then easily set a maximum 
wall temperature Tmax over which the efficient operation of a Hall thruster is not possible. 
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The temperature Tmax therefore represents a thermal constraint and it must be accounted 
for when designing a thruster. 

A semi-empirical time-dependent thermal model of a Hall effect thruster has recently been 
developed in order to determine the energy flux qp from a measurement of the temporal evo- 
lution of the channel wall temperature [16]. Yet, this model can be used the other way, i.e. 
to determine the wall temperature from the applied power and the channel sizes. Here, a 
simplified model of the thruster discharge chamber is used. The thermal enclosure is solely 
composed of the external and the internal cylindrical walls, as can be seen in Fig. [21 meaning 
that the anode and the rear part of the channel are not taken into account. A virtual cylinder 
that surrounds the enclosure is used to simulate the channel environment (coils, magnetic 
core...), as explained in [16]. In this work, only the steady-state wall temperature, i.e. the 
equilibrium temperature, is of relevance, meaning that the transient regime is ignored. More- 
over, our previous studies reveal that, to a great extent, heat conduction through walls can be 
neglected, only taking radiation heat transfer into account [16]. Under these assumptions, the 
total energy leaving a grey surface Ai is eiaTfAi, where a is the Stefan-Boltzmann constant, 
and the fraction of energy radiated towards another surface Aj can be calculated by multiply- 
ing the total radiated energy with the grey body configuration factor Gij [16]. The uniform 
temperature of the virtual cylinder is assumed to be T^nv = 450 K. As the temperature of the 
ceramic thruster walls is much higher, the virtual cylinder can be regarded as a mere heat 
sink. In order to assess the wall temperature, two quantities are needed: the energy flux qp 
as well as the thruster dimensions. The total power transferred to the channel walls Pioss was 
found to be of about 10% of the input power P for various Hall thruster [16]. It was also 
demonstrated that Pioss is almost split in equal amounts among the internal and external 
walls [TB] : 

(3.1) Ploss,int — AiYit X Qp,int — Pioss, ext — ^ext ^ Qp,ext — Pioss- 

Knowing the thruster dimensions L, h and d, the wall temperature T can be computed as a 
function of the input power. One must then check that the wall temperature does not ex- 
ceed Tmax- The validity of this simple thermal model has been put to the test by comparing 
the calculated and the measured wall temperature at a given applied power for several Hall 
thrusters. Results obtained with the PPSXOOO thruster are displayed in Table [TJ They clearly 
reveal a good agreement as the gap stays below 10 %. 

With BN-Si02 walls, a consistent value for T^ax is 970 K in compliance with outcomes of a 
study on the thermal behavior of Hall thrusters performed a few years ago |17j . The maximum 
authorized input power that corresponds to Tmax has been computed for a group of thrusters 
with different characteristic dimensions using a power loss factor of 10%. Results are shown 
in Table [2j Numbers are consistent with the range of operation power of tested thruster 
but for the micro-Hall and the NASA-built thruster. In the case of the micro-thruster, the 
computed maximum input power is only 9 W whereas this thruster operates at a power level of 
10 — 40 W [18]. This is not surprising as the thruster is equipped with a water cooling systems 
for the channel walls. The estimated maximum input power for the NASA-457M thruster 
is only of about 20 kW although this thruster is designed to be utilized around 50 kW. Yet, 
the origin of the discrepancy is well identified. In the thermal model, the power losses due 
to plasma wall interactions are supposed to be 10 % of the total input power for all thrusters 
without taking into account a possible size effect. As was shown in a previous work [1^, this 
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percentage is a function of the thruster size and it decreases when the size rises. In short, 
the trend originates from the fact that the ratio of surface (~ loss term) to volume (~ thrust 
generation term) decreases with the size. Indeed, when neglecting losses at the channel back, 
the surface to volume ratio is equal to for the channel of a HET. As a consequence, the 
percentage of power losses in the NASA-475M must be below 10 %. The thermal model gives 
a value of 5 % when the wall temperature Tmax is fixed to 970 K. 

Naturally, numbers given in Table [2] do not represent an upper limit in input power and a 
thruster can run at a power higher than the one given by the thermal model, as the BN-Si02 
ceramic can stand a temperatures larger than 970 K. The efficiency as well as the lifetime are 
nonetheless affected when the channel wall temperature is in excess of HOOK |20j . 

4 An empirical approach to the sizing method 

4.1 Description of the database 

A thorough open literature search using a wide range of resources combined with data- 
gathering within the French research program on electric propulsion allowed us to create 
a large database on Hall effect thrusters performances. The database contains information 
about thruster geometry as well as performances, notably the thrust T, the specific impulse 
Isp and the efficiency ry for a series of 23 different single-stage Hall thrusters. Moreover, the 
database includes information about the magnetic field strength 5, the discharge channel wall 
materials and the propellant gas. The entire database covers a vast range of input power that 
stretches from 10 W up to 50 kW and a large collection of data points in terms of applied 
discharge voltage and gas mass flow rate. A broad range of thrust level is certainly covered, 
going from 0.4 mN with a micro Hall thruster up to the 2.95 N delivers by the high-power 
thruster developed at NASA. The database also incorporates a few data about anode layer 
thrusters (TAL) of which the distinguishing features are a conducting channel wall and a 
short channel length [2]. 

A part of the collected data in terms of performance level is displayed in Fig. [3] and 
in Fig. m For all thrusters, channel walls are made of BN-Si02 and the propellant gas is 
xenon. The thrusters used to construct the two figures are the following: a 4 mm in diameter 
micro-Hall thruster operating at 10-40 W |18j . a laboratory model of the low power SPT20 
thruster [1^], a SPT50 thruster manufactured by the Kurchatov Institute [T^, the 1.5 kW- 
class PPS1350 HET developed and manufactured by Snecma [20], the 5kW-class PPSXOOO 
thruster which is a laboratory version of the PPS5000 technology demonstrator developed by 
Snecma [20], the 10 kW T220 designed and built by TRW and Space Power Inc. [21], as well as 
the 50 kW-class NASA-457M thruster [22]. The development of the thrust as a function of the 
discharge voltage is shown in Fig. [3] for the seven aforementioned HETs. The thrust of course 
increases with Ud- When operating at 0.2mg/s and Ud = 110 V, the micro-thruster delivers 
0.4 mN of thrust. On the opposite side of the thrust domain, the high-power NASA-457M 
thruster furnishes 970 mN of thrust when running at 35.2 mg/s and = 650 V. The evolution 
of the specific impulse along with the applied voltage is shown in Fig. HI The Igp increases 
with Ud and all thrusters follow an identical trend but the micro-thruster. The Snecma-built 
PPS1350 thruster delivers an Igp above 3250 s when it is fired at 1000 V in a low gas flow 
regime, as can be seen in Fig. [H An Igp of 3757 s was achieved by the SPT115 thruster at 
1110 V with a gas flow rate of 2.45 mg/s. The behavior of the anode thrust efficiency is shown 
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in Fig. m For most thrusters, a maximum is reached around Ud = 600 V when the walls are 
made of BN-Si02 ceramic. This specific behavior is likely to originate from the wall material 
properties combined with a change of the plasma properties at high voltage [20]. The plot in 
Fig. m reveals that the efficiency increases with the thruster size, as discussed in the preceding 
section. 



4.2 Validity of the scaling laws 

The propellant conversion efficiency a is the ratio of the ion mass flow rate to the propellant 
mass flow rate: 

4.1 a = — « — 

m Im 

Therefore a is not constant, as was considered in section 12. 2^ but it is a function of the 
discharge voltage and the mass flow rate. The value of a can be obtained from the performance 
data using the following equation: 



niiVi = amvi = am 



2eUd 



mi 



T I mi 

(4.2) ^a 



m y 2eU(i 

Figure [5] shows the calculated values of a for three different thrusters as a function of the 
applied voltage Ud when only Xe"*" ions are taken into account. This figure indicates that a 
depends both on the thruster size and on the value of Ud- The small SPT20 thruster indeed 
exhibits the lowest values of a. There is a natural link between r/, a and the dimensions. As 
can be seen in Fig. El ol increases quickly with the applied voltage, and for a large thruster it 
approaches unity at high voltage. The growth of a is especially connected with the electron 
temperature that increases with the voltage [23]. Actually, the production of multiply-charged 
ions must also be dealt for accurately assessing the value of a |24j . The growth of a is con- 
nected with both the the electron temperature and the fraction of multiply-charged ions [2l] . 
The two quantities increase with the voltage. According to the collected data set, the calcu- 
lated values of a vary between 0.3 and 0.96, not taking into account the micro thruster. For 
low voltages, a drops quickly due to a weak electron temperature. For an input power higher 
than 1 kW and an applied voltage above 300 V, the quantity a is commonly in the range 
0.8-0.9. The value determined by means of Eq. 14.21 are slightly underestimated as voltage 
losses are not accounted for. Nevertheless, the value of a obtained when using the ion velocity 
measured by way of a repulsing potential analyzer in the thruster near field are close to the 
ones computed with e ■ Ud kinetic energy [25]. That means the voltage losses term A is 
small. 

Figure [6] shows the thrust as the function of the product m^/Ud for a large collection of 
HETs and operating conditions. Apart from SPT20 thruster points at low power, all data 
points are aligned. As suggested earlier, the coefficient a is almost constant. From the slope 
of the curve in Fig. [6l one finds a = (0.89±0.01). Note that data about the SPT20 thruster 
are not considered for the linear fit. Equation 12.201 gives a relation between the thrust, the 
voltage and the characteristic dimensions. In Fig. [71 the thrust is plotted as a function of 
the product j^\/Udhd, as established by Eq. 12.201 for a variety of Hall thrusters firing with 
m ~ 5mg/s. The normal operating point of the thrusters ranges from 500 W to 5kW. As 
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can be seen, curves are linear, however, the slope depends upon the input power. A detailed 
analysis carried out with all thrusters available in the database reveals that thrusters group 
together depending on their normal operating power. Equation 12.201 was derived using the 
relation between m and the dimensions as well as the series of assumptions listed at the 
beginning of section 2.1. Though graph of Fig. [7] indicates that some of the assumptions may 
be too limiting. First, the propellant gas temperature is said to be constant an independent 
of Ud in such a way that the propellant speed inside the channel is fixed and it does not 
vary with Ud- In effect the gas temperature depends on the thruster geometry as well as 
on the input power. Second, the propellant conversion efficiency a, which was assumed to 
be constant, is a function of the discharge voltage as demonstrated hereinbefore. These two 
points are confirmed by the fact that the curves for thrusters of about the same sizes, e.g. 
the P5 and the PPSXOOO as weh as the PPS1350 and the SPTIOO, are close to each others. 
Yet, when the thrust T is plotted as a function of the product n„ hdy/Ud, all datapoints are 
aligned, as can be observed in Fig. [HI The plot is independent of the power level and the 
thruster sizes. This experimental fact reveals that the Melikov-Morozov criterion does not 
implies that the ratio \i/L is identical for all thrusters, whatever the geometry and design. 
It is even the opposite, and the mean value of the parameter (3 given in Table [3] must seen as 
the typical value that warrants a fine thruster functioning. 

As explained in section 12.11 the input power P naturally depends upon the characteristic 
thruster dimensions. Assuming Id ~ li, P is found to be a function of the ratio hd/L, see 
Eq. 12.211 However, using all gathered data, it appears that P varies linearly with the hd 
product, as shown in Fig. [9l This linear relation will be used in the remainder of the paper 
when sizing high-power Hall thrusters. 

At the end of section 12.11 two scaling laws have been established for the magnetic field 
B, see Eq. 12.131 and 12.151 The two plots in Fig. [10] show the evolution of B with l/L and 
1/hd, respectively, for five different Hall effect thrusters. As can be seen, data points are well 
aligned only when the hd product is considered. Therefore, one can consider that the scaling 
laws B oc 1/hd is more suitable than the one that solely incorporates the channel length L, 
although it was never considered up until then. 

Finally, a last remark is worth making. The atom density n„ that warrants an efficient 
functioning of a Hall effect thruster can be inferred from the xenon mass flow rate when 
the thruster dimensions are known, see Eq. 12.71 Assuming a propellant temperature of 800 K 
inside the channel and taking the gas mass flow rate at normal operating conditions, one finds 
n„ ~ 1.2 X 10^^ m~^ whatever the thruster picked in the database, but the micro-thruster. 
This number, which can be envisaged as an atom density constraint, is in good agreement 
with the value commonly found in literature. 

4.3 Determination of the proportionality coefficients 

In order to assess the required thruster dimensions by way of a scaling method for an available 
input power and a given thrust level, it is necessary to know the proportionality coefficients 
of equations encountered in section 2.2. As it was previously shown, the specific impulse is 
proportional to \/lJd- 




(4.3) 



with C/, 




= 123.4 



sp 
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This coefficient can also be inferred from the experimental function of \'T7d, see 

Fig.Sl The mean value of the coefficient determined by using the database is =99.7. 
The thrust is described by the following equation: 

T = CTi-rh y/Ud 

(4.4) with Ct, = aJ— = 1090.8, 

V rui 

when using a = 0.9. The coefficient Cti is the slope of the curve T versus m^/U2, see Fig. El 
Using the experimental dataset, one obtains = 1077.3. A second relation was established 
for the thrust: 

1 

I 




T = Ct2 -TVUdhd 

(4.5) with = OtninVnT^ 

Assuming that the electron temperature is 10 eV and the gas temperature is 800 K, the cor- 
responding velocities are f e = 2 x lO^m/s and Vn = 320 m/s, respectively. Furthermore, 
is taken to be 10^^ m~^ and ai is equal to 5 x 10"^'^ m^. Using those numbers, one finds 
= 7.65 X lO""^//?. The value of the (3 coefficient can be determined using our database as 
the size L is known for all thrusters. With numbers given above, one finds /3 = 0.007, and it 
comes = 0.109. The true experimental value C^^ is obtained by way of the type of curves 
plotted in Fig. [71 notwithstanding the fact that the slope of relation 12.201 depends upon the 
thruster normal operating power as discussed previously. A mean value is C^^ =0.0924. 

As shown in Fig. [9l the input power is a function of the product hd under our assumptions, 
and one can write: 

P = Cp- hd. 

(4.6) with Cp = ^^^«l.lxl0^ 

The value of Cp is calculated from the collection of operating points selecting various thrusters 
to vary the hd product. The empirical value of the proportionality coefficient is taken to be 
the slope of the line in Fig. [9l Cp = 1.2 x 10^. The value of all coefficients are summarized 
in Table O Finally, the thruster dimensions can be assessed afterwards, as it will be shown 
in the next section. 



5 Design of a high power Hall Effect Thruster 

High-power Hall effect thrusters in the range 10-30kW and able to deliver a thrust level 
around 1 N with a specific impulse of about 2500 s are thought to be used as the primary 
propulsion system for robotic space probes during interplanetary journeys [Tl 126^ [27]. Such 
high-power Hall thrusters may also be of interest for orbit transfer maneuvers of large satel- 
lites. Only a few high-power prototypes have been developed in the world so far and a sig- 
nificant research effort on this specific technology is now appearing within Europe. For these 
reasons, the sizing method based on aforementioned widely applicable scaling laws in combi- 
nation with our large database is employed to design a 20 kW-class thruster with T = 1 N as 
well as a 25kW-class Hall thruster with T = 1.5 N. 
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The discharge voltage of the two thrusters is fixed to Ud = 500 V, i.e. Isp = 2760 s when 
considering singly-charged ions and no voltage losses. Xenon is used as a propellant gas. The 
channel walls are assumed to be made of BN-Si02 ceramics. The highest wall temperature 
is set to 970 K. In this study, losses are fixed to 6 % of the applied power. The parameters a 
and /3 are set to 0.9 and 0.007, respectively, in agreement with previous analysis. 

5.1 Sizing of a 20kW-class thruster 

The process that consists in determining the thruster characteristic dimensions L, h and d as 
well as the magnetic field B must be carried out step by step. 

1. The required mass fiow rate is determined by means of Eq. 14.41 

2. The hd product is found using the relationship 14.61 

3. The channel length L is given by Eq 14.51 

4. The thermal constraint is used to assess the value of both the channel width h and the 
average channel diameter d. The process can be iterative, in the case the dimensions 
must be changed in order to satisfy the constraint: wall temperature below Tmax- 

5. The magnetic field strength is then obtained from the relation B oc 1/hd, see Fig. [TOl 

6. At last, it should be verified that the number density n„ is close to 1.2 x 10^^, as 
explained at the end of section 4.2. 

The channel length is the dimension that is the most difficult to determine as the proportion- 
ality coefficient Ct^ strongly depends on the value of f3. In fact, the ratio between Aj and L is 
not constant but it can vary considerably for different mass fiow rates, voltages and thruster 
types. Since the value of the hd product given by the method is quite reliable, it appears 
better to modify rather the channel length than the mean diameter in case an iterative loop 
is necessary to satisfy the thermal constraint. Indeed, outcomes of the thermal model are not 
much infiuenced by a change in the value of h as losses are fixed. 

For the 20 kW Hall thruster (IN, 500V), the sizing method leads to: m = 41.1 mg/s, 
L = 40.1mm, d = 250 mm and h = 66 mm and a magnetic field strength B = 136 G. The 
temperature is T^xt = 865 K and Tint = 962 K for the external, respectively the internal, 
dielectric channel wall. Finally, the value 1.12 x 10^^ m~^ is found for the atom number 
density. 

5.2 Sizing of a 25 kW-class thruster 

The same approach was applied to design a 25 kW Hall thruster (1.5 N, 500 V). The sizing 
method gives: m = 61.6 mg/s, L = 40 mm, d = 290 mm and h = 71mm. The magnetic 
field strength B is equal to 135 G. The temperature is T^xt = 871 K and Tint = 979 K for the 
external, respectively the internal, dielectric channel wall. The calculated number density is 
Un = 1.35 X lO^^m-3. 

The outcomes of our sizing method can be compared with results obtained in a study 
by M. Andrenucci and co-workers in which sizes are fixed instead of being computed [9]. 
They used their specific approach to design a 25 kW Hall thruster operating at 275 V. The 
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dimensions are fixed to L = 30 mm, d = 250 mm and h = 62 mm. They found a xenon mass 
flow rate of 88.7 mg/s, a thrust level of 1.49 N and a magnetic field of 135 G. Our procedure 
leads to a mass flow rate of 83.4 mg/s, an atom density of 2.61 x 10^^ m~^ and the same 
magnitude for B. Moreover, their thruster dimensions do not permit to respect our thermal 
constraint as T^xt = 970 K and Tint = 1103 K. Nevertheless, to achieve a thrust level of 1.49 N 
with Ud = 275 V, our approach predicts L = 35 mm, d = 310 mm and h = 67 mm, values that 
are not too far from the ones chosen a priori referring to a vast database [9]. 

6 Conclusion 

The Hall effect thruster sizing method described in these works considers the three charac- 
teristic thruster dimensions L, d and h, as well as the magnetic field strength B. The method 
relies on analytical laws that are established from the fundamental principles that govern 
the physics of a Hall thruster in the frame of simplifying assumptions. Besides, the thruster 
geometry must fulfill criteria about channel wall temperature and atomic number density. 
A vast database that incorporates 23 single-stage Hall thrusters covering a power range be- 
tween low and 50 kW allows to check the validity of scaling laws and to find the value of 
corresponding coefficients necessary to proportion a thruster. The sizing approach was then 
employed to obtain a proper estimate of characteristic dimensions and magnetic field strength 
for a 20 kW and a 25 kW Hall thruster capable of providing a thrust level of 1 N, respectively 
1.5 N. Our approach gave also satisfactory results when it was applied to check the design of 
a 100 W Hall thruster currently under development and test in our laboratory. 

Scaling laws developed here solely represent a first order approach due to the numerous 
simplifying assumptions on the physics at work in a Hall thruster. Nevertheless, for a given 
set of operating conditions they furnish a first estimate of the geometry and the magnetic 
field strength of a thruster, which permits to save time during the design and optimization 
stages. One way to improve the accuracy of the scaling method is to reduce the number 
of assumptions. The evolution of the electron temperature, of the gas temperature and of 
the fraction of multiply-charged ion species as a function of the discharge voltage could for 
instance be taken into account, however, available data are scarce that anyway limits the gain. 
Another way could consist in applying a statistical analysis of the vast database to directly 
obtain empirical sizing formulas. The design of experiment method could be an appropriate 
tool to achieve that goal. Finally, attempt to incorporate the magnetic field topology into 
scaling laws would represent a tremendous progress as the latter is the most fundamental 
feature to ensure a successful operation. 
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Table 1: Comparison between the measured PPSXOOO thruster external and internal wall 
temperatures and the ones obtained with a simple thermal model for three values of P. The 
level of power losses is fixed to 10%. 



Power (W) 


Text 


(K) 


Tint 


(K) 




Mes. 


Calc. 


Mes. 


Calc. 


1605 


690 


655 


740 


710 


2460 


762 


729 


788 


790 


4620 


833 


854 


862 


926 



Table 2: Outcomes of the thermal model in terms of maximum input power for various 
thruster types when T^ax is set to 900 K. The level of power losses is fixed to 10%. 



Thruster 


Maximum input 




power (W) 


micro 


9 


SPT20 


210 


SPT50 


870 


PPS1350 


2060 


SPTIOO 


2080 


PPSXOOO 


4200 


P5 


5550 


NASA-457M 


20300 



Table 3: Values of all proportionality coefficients needed for scaling laws. Values are in 
standard units. 



Coefficient 


Value 




123.4 


C*r 


99.7 


Cti 


1090.8 




1077.3 




0.109 




0.092 


Cp 


1.1 X 10^ 




1.2 X 10*^ 


p 


0.007 



19 



coils 



propellant 
supply 




cathode 



anode 



magnetic parts 



dielectric 
channel walls 



Figure 1: Schematic of a Hall effect thruster. The three characteristic dimensions L, h and d 
are also shown. 
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Figure 2: Drawing of the simple geometry used in the channel thermal model. 
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Figure 3: Thrust as a function of the discharge voltage for seven different Hall thrusters: 
micro-thruster (0.2 mg/s), SPT20 (0.472 mg/s), SPT50 (1.0 mg/s), PPS1350 (3.5 mg/s), 
PPSXOOO (5.0 mg/s), T220 (19.4 mg/s) and NASA-457M (35.2 mg/s). 
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Figure 4: Upper graph: Specific impulse as a function of the discharge voltage for seven 
different Hall thrusters. Lower graph: Anode thrust efficiency as a function of discharge 
voltage. The xenon mass flow rate is: micro-thruster (0.2 mg/s), SPT20 (0.472 mg/s), SPT50 
(1.0 mg/s), PPS1350 (3.5 mg/s), PPSXOOO (5.0 mg/s), T220 (19.4 mg/s) and NASA-457M 
(35.2 mg/s). 
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Figure 5: Propellant conversion efficiency a as a function of the discharge voltage for three 
types of Hall thrusters. The star symbol corresponds to values obtained for the PPSXOOO 
thruster when the ion velocity is recorded by means of a RPA. 
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Figure 6: Thrust as a function of rh^/Ud for a wide range of thruster sizes. 
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Figure 7: Thrust as a function of the product j^\/Udhd for the SPT80 at m = 4.58 mg/s as 
well as the SPTIOO and the P5 at ttt. = 5 mg/s. 
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Figure 8: Thrust as a function of nny/Udhd for different thrusters; Upper Graph: SPT80 
(4.58 mg/s), SPTIOO, PPS1350i (4.82 mg/s), PPS13502 (5.21 mg/s), PPSXOOO and P5i at a 
flow rate of about 5 mg/s; Lower Graph: SPT80 (4.58 mg/s), SPTIOO (5 mg/s), PPS1350i 
(4.82 mg/s), PPS13502 (5.21 mg/s), PPSXOOO (5 mg/s), P5i (5 mg/s), P52 (10 mg/s), T220i 
(15.9mg/s), T2202 (19.4mg/s), NASA-457Mi (20.0mg/s) and NASA-457M2 (24.9mg/s). 
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Figure 9: Input power as a function of the product h x d for six Hall thrusters. Also shown 
is a linear fit through all data points. 
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Figure 10: Upper graph: Magnetic field strength as a function of 1/L. Lower graph: Magnetic 
field strength as a function of 1/hd. 
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